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We investigate the one-loop supersymmetric QCD (SUSY-QCD) and electroweak (SUSY-EW) 
corrections to the top quark decay into a 6-quark and a longitudinal or transverse IF-boson. The 
corrections are presented in terms of the longitudinal ratio T(t — > Wz,6)/r(t — > Wb) and the trans- 
verse ratio F(t — > W-b)/T(t — > Wb). In most of the parameter space, both SUSY-QCD and 
SUSY-EW corrections to these ratios are found to be less than 1% in magnitude and they tend to 
have opposite signs. The corrections to the total width F(t — > Wb) are also presented for comparison 
' with the existing results in the literature. We find that our SUSY-EW corrections to the total width 

^vQ , differ significantly from previous studies: the previous studies give a large correction of more than 

10% in magnitude for a large part of the parameter space while our results reach only few percent 
at most. 
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■ I. INTRODUCTION 

00 ' 

As the most massive fermion in the Standard Model (SM) , the top quark is believed to be sensitive to new physics 
[1]. Since the measurements of top quark properties in Run 1 at the Fermilab Tevatron have only small statistics, 
there remains plenty of room for new physics to be discovered in future precision measurements. Such possibilities 
exist on our horizon: while the upgraded Tevatron collider will significantly improve the precision of the top quark 
measurements, the CERN Large Hadron Collider (LHC) and the planned Next-generation Linear Collider (NLC) will 
serve as top quark factories and allow precision measurements of the top quark properties [2] . 

Although various exotic production processes [3] and rare decay modes [4] of the top quark will serve as robust 
evidence for new physics (since these exotic processes are forbidden or severely suppressed in the SM), the role of 
the dominant decay mode t — » Wb in probing new physics should not be underestimated [5] . One advantage of this 
decay mode is that it is free of non-perturbative theoretical uncertainties [6] and future precision experimental data 
can therefore be compared with accurate theoretical predictions. The other advantage of this channel is that the 
VF-boson as decay product is strongly polarized and the helicity contents (transverse-plus W+, transverse- minus W— 
and longitudinal Wl) of the VF-boson can be probed through a measurement of the shape of the lepton spectrum in 
the TF-boson decay [7]. Therefore, the study of top decay into a polarized FF-boson will provide additional information 
about the tWb coupling 1 and has so far attracted a lot of attention from both theorists and experimentalists. On 
the experimental side, the CDF collaboration has already performed the measurement of the helicity component of 
the IF-boson in top quark decay from Run 1 data and obtained the results 

r L /r = 0.91 ± 0.37(stat.) ± 0A3(syst.), 

r+/r = o.ii ±0.15 , 

where T is the total decay rate of t — > Wb, and Tl and T + denote respectively the rates into a longitudinal and 
transverse-plus VF-boson. Although the error of these measurements is quite large at the present time, it is expected 
to be reduced significantly during Run 2 of the Tevatron and may reach 1% ~ 2% at the LHC [8]. On the theoretical 
side, predictions of these quantities in the SM up to one-loop level are now available [9,10]. For the mass values of 
Refs. [9,10], the tree-level results are 0.703 for T L /T, 0.297 for T^/T and e>(10~ 4 ) for T + /r, and the QCD corrections 



1 Among the three polarizations of the IF-boson in top decay, the longitudinal mode is of particular interest since it is relevant 
for the understanding of the electroweak symmetry breaking mechanism. 
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to these predictions are respectively —1.07%, 2.19% and 0.10%, while the electroweak corrections are at the level of 
a few per mill. 

In order to probe new physics from the future precise measurement of Tl/T, T_/r or r + /T, we must know 
the new physics contributions to these quantities in various models. In this article we focus on the supersymmetric 
contributions in the framework of the Minimal Supersymmetric Model (MSSM). In this model the one-loop corrections 
to the total width of t — > bW have been studied by several groups [11-15]. However, the corrections to Tl/T, T_/r 
or r+ /r arc still lacking and calculating these corrections is the major goal of this article. 

In the MSSM the genuine supersymmetric corrections (i.e., the corrections from various sparticle loops) are SUSY 
electroweak (SUSY-EW) corrections, arising from interactions of charginos or neutralinos, and SUSY-QCD corrections, 
arising from interactions of gluinos 2 . The SUSY-EW and SUSY-QCD corrections to the total width were first 
investigated in [11,12] without considering squark mixings, and later were re-studied extensively in [13-15] by taking 
into account the squark mixings. Although our aim in this paper is to study the corrections to Tl/T, T-/T and 
r + /r, we will also repeat the calculations of the corrections to the total width for two reasons: One of the reasons 
is that the previous SUSY-EW calculations were performed many years ago and some SUSY parameter space viable 
at that time has been ruled out by recent experiments. The other reason is that the SUSY-EW corrections to the 
total width were found to be exceptionally large ( exceeds 10% in a large part of parameter space) [13] and should, 
therefore, be checked carefully. In fact we find they are much smaller for reasons we will discuss. 

This paper is organized as follows. In section II, we present relevant formula for our calculations. In section III, 
SUSY-EW corrections are calculated and compared with previous studies. In section IV, we examine SUSY-QCD 
corrections. For the numerical calculations of both SUSY-EW and SUSY-QCD corrections, we first perform a scan 
in the typical allowed parameter space to find out the typical size of the corrections, and then consider some special 
scenarios (such as a very light sbottom or gluino) to figure out if the corrections are exceptionally large in some regions 
of parameter space. The conclusions are given in section V. 



II. FORMALISM 



In this section we present the formulas for calculating new physics contributions to t — ► bW in the so-called "G_f- 
scheme" [17] where the Fermi constant Gf and the pole masses mw, mt, m b, • • • arc chosen as input parameters. This 
on-shell scheme is the most convenient one for studying such new physics contributions. These formulas are valid for 
studying one-loop corrections to the decay t — > Wb in any renormalizable new physics model. 

The tWb vertex at one-loop level takes the form [17] 

r p = -z2 3 ^Glm w V tb ^P L (l + ±5Zt + ^6Z t L + 5Z? - 8Zf - l -U w - ±Ar™ + F L ) 

+1^PrFr + p^P L H L + p^PrHr} . (1) 

Here Pr,l = \ (1 ± 7s) are the chirality projectors. The form factors Fl,r and Hl.r represent the contributions from 
the irreducible vertex loops. (For the comparison of our results with Ref. [13], we adopted the same notations for the 
form factors as in Rcf. [13].) SZ^ and 5Z^ denote respectively the field renormalization constant for bottom quark 
and top quark. SZ^ 2 are the renormalization constants for the SU(2) gauge fields W® and the coupling constant g 2 , 
defined as W° -» {Z¥) 1/2 W« and g 2 -» Z^ (Z 2 w )- 3 / 2 g 2 with Zf <2 = 1 + bZf 2 . Their difference is given by [17] 

5ZY - 5Zf = — ^2^1 (2 ) 

swcw mz 

where Cw = cos^w, sw = sin^v^i and S 7 ^(0) is the unrenormalized 7-Z transfer at zero momentum. Tl\y in Eq.(l) 
represents the finite wave function renormalization of the external lU-boson and is given by = E^(m^) + SZ^ , 
where Y,' w = dT,w(p 2 )/dp 2 with denoting the VF-boson self-energy. Ar now is the corresponding new physics 
contribution to Ar, which is the radiative correction to the relation of Mw, Mz, a (the fine structure constant) and 
G F defined by [17] 



2 The corrections from the Higgs interactions to the total width of t — » bW are not unique to the MSSM since they also exist 
in the general two-Higgs-doublet models. Such corrections are found to be less than 1% for the whole parameter space [16]. 
Therefore we do not calculate such corrections in this paper. 
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By decomposing Ar now into the universal part and non- universal part, i.e., Ar ncw = Arg° w + Ar^, and inserting 
the explicit forms of 5Z^ 2 , one can re-express T M as [17] 



-hl' w {m 2 w ) - l- Ew{0) 2 ^ - ^Ar^ + *l) + 7^*fl +Pt„P L H L + p^P R H R ) . (4) 

In our calculations we adopted all the conventions in [17] except for the fermion-field renormalization. To be valid in 
general cases, the self-energy of a fermion is decomposed as 

S/(p) = ^(p 2 ) iP L + Ef (p 2 ) iP R + m f Y,f{p 2 )P L + m f X s f R (p 2 )P R . (5) 

Using the on-shell conditions for external fermion fields [18], we obtain SZf as 



SZf = Zj(mj) + mj [Ef(ro/) + S / ( m f) + S / K) + £/(«*/) - ( 6 ) 
<SZf = (m 2 f ) + m 2 \Z?(m 2 f ) + tf(m 2 f ) + ^ L '{m 2 ) + tf R '(m 2 )\ , (7) 

where for S one takes the real part of the loop integrals only but keeping all parts for the coupling constants in the 
fermion self-energies. 

The rate of the top quark decay into a polarized TU-boson can be obtained by using the explicit expression of 
polarization vector of IU-boson 3 or using the project technique introduced in [9,10] . Through some tedious calculations 
we obtain 

r L = G F m w™t\v tb \ 2 (i-x 2 ) 2 ^ + Re{sz ^ + + 2gz w _ 2Sz w _ Uw _ Ar , ncw + 2Fl) 

8v 2ir x 

+Re(H R )m t (l - x 2 )} , (8) 

T_ = G F m wJh\VtA 2 (1 _ 3.2^2^ + Re{SZ L + 5Z L + 2SZ W _ 2§Z W _ ^ _ Af » + ^ (Q) 

where (r_) denotes the rate of the top quark decay into longitudinal (transverse-minus) IF-boson and x = mw/m t . 
In deriving Eqs.(8,9) we have neglected the 6-quark mass for simplicity and this will introduce an uncertainty of several 
permille on Tl.-- Another consequence of neglecting m b is T + = due to angular momentum conservation [9] and, 
as a result, the total decay rate of t — > bW is obtained by Y = Yl + r_. In our following calculations, we retain the 
bottom quark mass only when it appears in couplings or in the sbottom mass matrix [19] since in those cases the 
bottom quark mass is multiplied by tan/? and can not be neglected for large tan/3. 
We define the ratios 

f il _ = r ii _/r ) (io) 

which can be measured in experiments. In our numerical results we will present the relative SUSY corrections to 
them, i.e., 5Tl_/Y° l _ with STl,- denoting the SUSY corrections and r° denoting the SM predictions. 

III. SUSY ELECTROWEAK CORRECTIONS 

In this section we investigate the SUSY clectroweak corrections. First we consider vertex corrections and quark 
self-energy corrections as depicted in Fig. 1. This part of the corrections are expected to be sizable for large tan j3. In 
Appendix A, we list in detail the relevant Feynman rules and present our analytical results. We checked that our results 



3 For the construction of covariant polarization vector of IF-boson, see Eq.(5.16) in [18]. 
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do reduce to those in [11] when switching off squark mixings and neglecting the &-quark mass in quark-squark-chargino 
(or neutralino) interactions. 

The SUSY-EW contributions to the self-energies of gauge bosons, which are necessary to calculate for SZ^ 2 and 
Hw, have been calculated by several groups [11,20]. We recalculated them and find that our results agree with those 
in [11] when switching off squark mixings. In our numerical calculations we will take into account the mixings between 
top-squarks and between sbottoms. 

As to the SUSY-EW contributions to Ar, both the universal and non-universal parts are nicely presented in [21]. 
We checked the analytic expressions and incorporated them in our FORTRAN code. Our numerical results show that 
the total contribution to the width T(t — > Wb) from 5Z^ 2l U w and Ar SUSY-EW is generally less than 1%. 




bs (t a 
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FIG. 1. Some of the Feynman diagrams of SUSY-EW corrections to t — > Wb: (a-c) are vertex loops; (d,e) are top and bottom 
quark self-energy loops. 
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FIG. 2. SUSY electroweak corrections to the width T(t — > Wb) versus sbottom mass. The solid curve is our result and the 
dashed curve is taken from Fig. 7 of Ref. [13] . The input parameters are those of case C in Fig. 7 in Ref. [13] . 



As pointed out in Appendix A, our vertex corrections are quite different from those in [13]. To see the numerical 
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difference between our resuits and those in [13], we present a comparison in Fig. 2, which corresponds to the case C 
of Fig. 7 in [13]. We see from Fig. 2 that the difference is significant. The results of [13] are much larger in magnitude 
than ours: in a large part of parameter space the results of [13] exceed 10% in size while our results are below 5%. 

Fig. 2 also shows a common feature for both results, namely, there exists a peak in the correction at m~ h ~ 109 
GeV. This peak comes from the "threshold effect" [22], which is explained in Appendix B. 

Note that some of the input parameters in Fig. 2, which were allowed at the time when the calculations in [13] were 
performed, may no longer be allowed nowadays. In our following numerical calculations, we will consider the current 
constraints on the parameters. 

The relevant SUSY parameters in our calculations are: (1) tan/3, Mq, Mjj, Md, P and A t .b which enter the mass 
matrices of top-squark and sbottom [19] (see Appendix A); (2) soft-breaking gaugino masses M\ and M 2 appearing 
in the chargino and neutralino mass matrices (see Appendix A) ; (3) the masses of sleptons and the squarks of the first 
two generations involved in the calculation of SZ^ 2 , and Ar SUSY_EW . To lessen the number of the parameters in 

(3) , we assume an universal soft breaking parameter Ml and neglect left-right mixings for sleptons and the squarks of 
the first two generations. Such assumptions do not affect our numerical results significantly since, as we said before, 
the contributions from SZ^ 2 , Tlw and Ar new are small for the allowed parameter space. As for the soft breaking 
gaugino mass parameters M\ and M2, we assume the relation Mi ~ Q.5M2 [23]. With these assumptions the relevant 
SUSY parameters are reduced to 

fhQ,mu,friD,At,Ab,fi,ta,ni3,ML,M2. (11) 

To estimate the size of the SUSY-EW effects, we first performed a scan over the nine parameters in Eq. (11). 
In our scan we restricted the parameters with mass dimensions to be less than 1 TeV and considered the following 
experimental constraints: 

(1) \i > and a large tan/3 in the range 5 < tan/3 < 70, which seems to be favored by the muon g — 2 measure- 
ment [24]. 

(2) 5p < 0.002 [25], which will constrain the mass splitting between stops and sbottoms. (We use the program 
FeynHiggs [26] to generate the values of 8p and the Higgs masses.) 

(3) The LEP and CDF lower mass bounds on sparticles and the lightest CP-even Higgs boson [25,27] 

mi > 86.4 GeV, m~ b > 75.0 GeV, > 67.7 GeV, m^o > 40 GeV, 

mj > 95 GeV, mp > 41 GeV, m~ q > 138 GeV, m h o > 91 GeV, (12) 

where rriq denotes the mass of squarks of the first two generations. 
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FIG. 3. The scattered plot in the plane of ST/T versus ST L /r L . 
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Our scan results are shown in Fig. 3 in the plane of <5r/r° versus 8Tl/T q l , In our calculations throughout this paper, 
we have fixed G F = 1.16637 x 1(T 5 GeV~ 2 , m w = 80.451 GeV, m t = 174.3 GeV and m b = 4.8 GeV. Fig. 3 shows 
that the SUSY-EW correction is generally below one percent in size and tends to be positive for T and negative for 

Although the SUSY-EW correction turns out to be quite small (below one percent in size) in the dominant part 
of the allowed parameter space, we found from our scan that the correction can reach a few percent in some narrow 
corners of the parameter space. These corners have two feathers: (1) /j, <C M2 so that the lighter chargino and 
the lightest neutralino are Higgsino-like, and tan (3 is large so that the 6-quark Yukawa couplings in quark-squark- 
Higgsino interactions can be greatly enhanced [19]; (2) The two sparticles in a self-energy loop of top quark lie below 
the threshold, or in other words, the sum of the masses of the two sparticles is lighter than top quark mass. Generally 
speaking, for m fc ~ + ni_^+ < m t the sbottom^! )-chargino(x^~) loop can give a relatively large contribution. With the 
increase of m fe ~ and , the contributions from this loop decrease as required by a decoupling theorem. A peak in 
correction size occurs at rog i + — m tj &s shown in Fig. 2 and explained in Appendix B. To show the size of the 
corrections in such corners of the parameter space, we present some numerical results in two special scenarios: 



Scenario I: We assume \x <C M 2 and assume an universal soft-breaking mass M 



SUSY 



for squarks. 



• Scenario II: We assume a very light sbottom b\ of about 5 GeV. So far such a light sbottom has not been ruled 
out by current experiments if the sbottom mixing angle 0& is tuned to satisfy cos6>(, ~ 0.38 [28]. For such a light 
sbottom, tan (3 must be large enough to cause a sufficiently large mass splitting between the two sbottom mass 
eigenstates. In addition, the lighter top-squark t\ cannot be heavier than 300 GeV in order to satisfy the 5p 
bound. 

The results in Scenario I are plotted in Figs. 4 and 5. We present our result only in the region where the constraints 
given in the paragraph following Eq.(ll) are satisfied. The peaks in these figures happen at + = mt where 
and m~+ are the lighter sbottom and top-squark masses, respectively. Fig. 4 shows that below the threshold 
the SUSY-EW corrections to the width can be as large as 3% for tan/3 > 50 and after crossing the threshold, the 
corrections change sharply to —8% and then quickly decrease in size. Concerning the behavior at threshold, we want 
to point out that the results within the region < + m% — m t < T/2 are not reliable since perturbative expansion 
of the ^-matrix element breaks down for these points [13] (also see discussions in Appendix B). Fig. 5 shows that 
<5r_/r°_ can reach 1% near the threshold and its size is generally larger than that of STl/T 1 ^ . 
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FIG. 4. SUSY electroweak corrections to V versus Msusy in Scenario I. 
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FIG. 5. SUSY eletroweak corrections to T-,l and F_/ L (= f_/fi) versus Msusy in Scenario I. 

Now we turn to Scenario II. For numerical calculations in Scenario II, it is inconvenient to take the parameters in 
Eq.(ll) as input. Instead, we may choose 



m~ bl , m fi , mj 2 , (i, A b , 9 b , tan/3, M L ,M 2 



(13) 



as input parameters. Using the formula listed in Appendix A, one may relate these two sets of parameters. In our 
numerical calculations, we keep /z and tan /3 as variables and study the dependence on them. Other parameters are 
fixed to be 4 



m h = 150 GeV, m h = 500 GeV, M L = A b = 300 GeV, 
m~ bi = 5 GeV, cos(9 h = 0.38, M 2 = 400 GeV. 



(14) 



Note that we fixed such values for mg i and cos b since this scenario is characterized by assuming a very light sbottom 

b\ of about 5 GeV with the sbottom mixing angle sytisfying cos#b ~ 0.38 (see the definition of this scenario). 

The results in Scenario II are plotted in Figs. 6 and 7. In Fig. 6 we plot ST /T° versus /i and in Fig. 7 we plot 
<5r_,L/r°_ i and 5r_/ L /T a _^ L versus [i. Again we see that a large tan/3 can cause large SUSY-EW corrections. Fig. 6 
shows a broad region of fi in which corrections to the width can be a few percent. Like in Fig. 4 the peaks occur at 
m bi + rn~T = m t . Fig. 7 shows that, like in Scenario I, the SUSY-EW corrections to T_ and Tl are of opposite sign 
and the size of the former is larger than that of the latter. 



4 We found from numerical calculations that the results are not sensitive to the parameters m^, wit 2 , Ml and Ab. 
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IV. SUSY-QCD CORRECTIONS 



In this section we investigate SUSY-QCD effects. Compared with SUSY-EW corrections, SUSY-QCD corrections 
are easier to calculate since SZ^ 2 , and <5r new receive no contribution from SUSY-QCD interaction at the one- loop 
level and consequently, one only needs to calculate quark self-energy and vertex correction depicted in Fig. 8. 




FIG. 8. Feynman diagrams for SUSY-QCD contribution to t — > Wb. 

Our analytical results for the form factors -Ft.fi and Hl.r agree with those in Ref. [14] and thus we do not present 
our analytic results here. Our numerical results for the corrections to the total width T(t — > Wb) agree with Ref. [15] 
for the same choice of parameters. 

The parameters involved in numerical evaluations are the gluino mass nig, and the parameters toq, %, rhu, A t , Ab, 
fi and tan f3 which enter the mass matrices of top-squark and sbottom [19] . In order to compare the size of SUSY-QCD 
corrections with that of SUSY-EW corrections, we have shown our scan results of SUSY-QCD corrections already 
in Fig. 3 in the preceding section. In our scan we have considered the relevant constraints given in the paragraph 
following Eq.(ll) and we also required nig > 200 GeV, which is the current experimental bound for the gluino in 
the framework of the minimal supergravity model [25]. From Fig. 3 one sees that the SUSY-QCD corrections and 
SUSY-EW corrections tend to have opposite signs and thus may cancel in large parts of parameter space. This differs 
drastically from the conclusion of Ref. [13,14] which claims that the two type corrections have the same sign and thus 
the combined contribution can reduce the width T(t — > Wb) by as much as 25%. 

We see from Fig. 3 that in the allowed parameter space with nig > 200 GeV the SUSY-QCD corrections are smaller 
than 1% in magnitude. However, we found that the corrections can reach a few percent in some corners of the 
parameter space if we relax the bound nig > 200 GeV (this bound is model-dependent). In the following we present 
some numerical results by relaxing this bound nig > 200 GeV, and in order to compare with SUSY-EW corrections 
we consider the same two scenarios as in the preceding section. 

The results in Scenario I are plotted in Figs. 9 and 10 for different m§ values in the allowed region of Msusy- We 
fixed u, — 80 GeV and tan j3 = 70, and included the possibility of a very light gluino of 16 GeV, which has not been 
excluded by current experiments 5 (for comparison we also plot a curve with a 100 GeV gluino although it may be 
severely disfavored by existing experiments). Fig. 9 shows that the magnitude of SUSY-QCD corrections decrease as 
Msusy or nig becomes large. Comparing with the SUSY-EW corrections in Fig. 4, one sees that the maximum size 
of the SUSY-QCD corrections is smaller than that of the SUSY-EW corrections, which can be greatly enhanced by 
the 6-quark Yukawa coupling for large tan (3. 



5 Such a very light gluino has not been excluded by experiments and was recently used to solve the long-standing fo-quark 
production puzzle [29]. To achieve that, a light sbottom of about 5 GeV is also required, which leads to a stronger constraint 
from Rb on the SUSY parameters [30] . 
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FIG. 9. SUSY-QCD correction to T versus Msusy in Scenario I. 
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FIG. 10. SUSY-QCD correction to F-,l and F_/l versus Msusy in Scenario I. 
Next we present results in Scenario II. Similar to the SUSY-EW case, in Scenario II it is convenient to choose 

TO &! ' m h > m i 2 ) Mi #6) tan P, m g (15) 



10 



as input parameters. Our numerical studies show that the SUSY-QCD corrections are only sensitive to m^ and rrig. 
So we fix m~ bi = 5 GeV, m t ~ a = 500 GeV, n = 80 GeV, A b = 300 GeV, tan/3 = 70 and cos6> 6 = 0.38 and plot the 
corrections in Figs. 11 and 12 for different gluino masses. Since in this scenario we assume a very light sbottom 
(m^ = 5 GeV), we do not consider the possibility of a very light gluino (m g — 16 GeV) due to the Rb constraints [30]. 
From the figures we see that for a m^ — 100 GeV and m~ g — 100 GeV, SUSY-QCD corrections can reach —3%, 
comparable in magnitude with SUSY-EW corrections. 




M~ t (GeV) 

FIG. 11. SUSY-QCD corrections to T verus the lighter top-squark mass in Scenario II. 
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FIG. 12. SUSY-QCD corrections to F-,l and F_/l versus the lighter top-squark mass in Scenario II. 
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V. SUMMARY AND CONCLUSION 



We investigated the one-loop SUSY-QCD and SUSY-EW corrections to the top quark decay into a 6-quark 
and a longitudinal or transverse lU-boson. The corrections are presented in terms of the longitudinal ratio 
T(t — ► WLb)/T(t — ► Wb) and the transverse ratio T(t — ► W_b)/T(i — » Wi>). In order to compare with the exist- 
ing results in the literature, we also present the corrections to the total width T(t — > W£>). To find out the typical size 
of these corrections, we performed a scan in the typical allowed parameter space. In order to decide if the corrections 
are exceptionally large in certain corners of the parameter space, we then considered some special scenarios (such as 
a very light sbottom or gluino). Our observations are: 

(1) In most of the parameter space, both SUSY-QCD and SUSY-EW corrections to the ratios are less than 1% in 
magnitude and they tend to have opposite signs. Only in some corners of the parameter space the corrections 
can reach a few percent. 

(2) The corrections to the total width T(t — ► Wb) agree well with previous calculations for the SUSY-QCD part, 
but differ significantly from previous calculations for the SUSY-EW part. Unlike the previous studies, which 
showed a large SUSY-EW correction of more than 10% in magnitude for a large part of the parameter space, 
our SUSY-EW results are only a few percent at most. 

We conclude by remarking that these SUSY corrections, despite their smallness in size, should be taken into account 
in the future precision tests of the top quark decay t — > Wb to determine whether SUSY is indeed the new physics 
chosen by Nature. (As far as we know, the corrections to the tWb coupling in other kinds of new physics models 
are also very small [31].) However, from the viewpoint of probing SUSY through revealing its effects in the coupling 
tWb, the decay t — > Wb may not be the best channel. Some single top quark production channels proceeding through 
the tWb coupling may be complementary or even do better in this aspect. For example, SUSY effects in the tWb 
coupling may cause observable effects in single top quark production in hadron collisions [32], and the single top 
quark production via electron-photon collision at a future linear collider is also quite sensitive to anomalous tWb 
coupling [33]. 
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APPENDIX A: SUSY-EW LOOP RESULTS 

In this appendix we first list the mass matrices of sfermions, charginos and neutralinos involved in our calculations, 
then write down the relevant interaction Lagrangian and finally present the explicit expressions of the SUSY-EW 
contributions from the loop diagrams shown in Fig.l. 

1. Mass matrices of sfermions, charginos and neutralinos 

Assuming no generation mixing for squarks in the soft-breaking terms, the mass-square matrices for up-type and 
down-type squarks in each generation are given respectively by [19] 

M 2 = ( + m z cos2 P(h ~ eusm 2 e w ) + to 2 m u (A u - fj, cot (3) 

u \ m u (A u - fj, cot f3) Mjj + m| cos 2(3e u sin 2 6 W + m 2 u 

M 2 _ ( m q ~ m|cos2/3(i + e d sin 2 9 W ) + m 2 , m d (A d - ,utan/3) 

<* \ m d (A d - ytitan/3) M|> + m| cos 2(3e d sin 2 9 W + m 2 , 

where Mq, Mu and are soft-breaking mass terms for left-handed squark doublet, right-handed up and down 
squarks, respectively. A u (A d ) is the coefficient of the trilinear term H 2 QU (HiQD) in soft-breaking terms, and 
tan/3 = v 2 /vi is the ratio of the vacuum expectation values of the two Higgs doublets. 
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For the first two generation squarks, the left-right mixings can be neglected due to the smallness of relevant quark 
masses and thus the weak-eigenstates are mass-eigenstates. From Eqs.(Al,A2), we then get the relation between the 
left-handed up-type squark mass M^ L and the left-handed down-type squark mass Mj l 

M? L = Mj L + m| cos 2/3(1 - sin 2 9 W ) + m 2 - to 2 , . (A3) 

For the third generation squarks, the left-right mixings cannot be neglected and the mass-eigenstates q lt2 (q = t, b) 
are related to the weak-eigenstates q_h,R by an unitary rotation R 

1 1 )=r(I l ) = ( COsd l Sin ^ )(l L ) , (A4) 
qi ) \qn ) ^-sin0q cos9 q J \q R J 

where 9 q is the so-called mixing angle between and q~ R . 

The mass-square matrices in Eqs.(Al, A2) can be alternatively expressed by the squark masses and mixing angle 

as 

M 2 = ( C ° s2 dqm ^ + Sin2 *« m fc S[n0q cosd i( m i ~ m V \ ( A5 ) 
q \ s\n9 q cos9 q {m 2 qi - m|J sin 2 9 q m\ x + cos 2 9 q m\^ ) 

By comparing Eq.(A5) with Eqs.(Al, A2), one can obtain the relationship between the two set of input parameters 
in the squark sector, {m qil m q21 9 q ) and (Mq, Mjj, ^d, M.b, m)- The mass matrices for sleptons take the similar 
forms as for squarks. But due to the smallness of lepton masses, the left-right mixings can be neglected for the three 
generations. 

The mass matrices of charginos and neutralinos are also involved in our calculations. The chargino mass matrix 
can be diagonalized by two unitary matrices U and V 



U* 



f- )V = diag^TO~+f (i = l,2). A6 

miyv2cosp jj, J I *< J 



The neutralino mass matrix is given by 

/ Mi — mz sin9w cos/3 mz sin 9y/ sin /3 \ 

Y_ Mi mz cos9w cos f3 — to^ cos #w sin/3 

— mz sin 9w cos/3 mz cos 9\y cos /3 —fi 

\ m z sin 9 W sin — m z cos 9 W sin /3 — /j, / 



(A7) 



and it can be diagonalized by an unitary matrix N 

^YN- 1 = diagjm^o} {j = 1,2,3,4). (A8) 

In the above expressions Mi and M 2 are respectively the soft-breaking U(l) and SU{2) gaugino masses. m~+ and 
m^o are the masses of charginos and neutralinos, respectively. One point we want to mention is that U, V and N 
defined in Eq.(A6, A8) are not unique and one should choose U, V and A in a way that all the masses of charginos 
and neutralinos are positive. 

2. Interaction Lagrangian 

The interaction Lagrangian relevant for our calculations is given by 

C = -f {t* a X7«i B* ail5 )b + bfxt(A% - B b 0i75 )t + t*X{A«f - Blf^t + bfx^f - B h ^ l5 )b] 
-^=i£ii$i** & b f3 W+ + gfj^O^PL + 0*P R )x-W+ 

+y + Rki^Jk^W-W^ + h.c. , (A9) 
where the sum over repeated indices is implied, and 
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AL = RUV^ - X b U t2 ) - R l a2 \ t v: 2 , (A10) 

Bli = KiiVa + X b U l2 ) - R^XtV* , (All) 

A b (jl = i?^(C/* - \ t V l2 ) - R h f32 \ b U* 2 , (A12) 

B% = RU U n + X tV l2 ) - R b 02 X b U* 2 , (A13) 

= Kii^Nji + V^^N*, tandw + ±N* 2 )) + < 2 (A t JV* 4 - V^-N^ tan%) , (A14) 

<° } = RU-*tN j4 + yft^Nfr tzn8 w + ±N* 2 j) + RU^tN* 4 + y/2^N dl tan%) , (A15) 

1 1 p2i 

A T = R b pi (X b N j3 + V^i-N*, tan0 w - -N* 2 )) + R b p2 (X b N* 3 + ^-N }1 t&n0 w ) , (A16) 

KT = R (}i(- X b N 33 + ^(\ N Ji ^9 W - ±N* 2 )) + R b p 2 (\ b N* 3 - ^-N n tan9 w ) , (A17) 

0% = ~^ffM + ^Wi , (A18) 

0^ = ^=N* 3 U i2 + N* 2 Ua . (A19) 

Here X t = mt . . and At = 7= " lb -. As for the notations for chargino and ncutralino mass matrices, we follow 



the notations of Ref. [19] rather than those of Ref. [13]. Note, however, that the resulting interactions are the same. 

3. Analytic results of SUSY-EW loops 

To compare our results with those in Ref. [13], we present our results in the same way as in Ref. [13]. First, from 
matrices in Eqs.(A10-A19) we define the following matrices 

A t ± =A t ±B t , Af ] = A^ ± B*(°) , (A20) 
A b ± = A b ± B b , Af } = A b ^ ± B b ^ , (A21) 

and construct the combinations as 



A (l) 


= -A^C^A*™ 


£7(1) 


= -A t *O R A t l 0) , 


(A22) 




= -A^A^ 




= -A t *0 R Af ) , 


(A23) 


Cr(i) 


= -A l ;o L A m 


G (D 


= -A u O L A m , 


(A24) 


D (i) 


= -A l ;o L Af ] 




= -A t *0 L Af ) , 


(A25) 


A (2) 


- fi>' U b A m 


* A t(o) 


C'V = -R^R^AV* A m , 


(A26) 


£>(2) 


T>t r>b MO) 

— n al itp 1 ^i + 


* aKo) 
A + > 


= -R^R^A^A^ , 


(A27) 


A (3) 


= Af^O^t , 


-- 


= A m *0 L A b _ , 


(A28) 


B (3) 


= Af h O L A b + , 


F (3) 


= A^'C^A*. , 


(A29) 


C (3) 


= Af ) *0 R A b _ , 


-- 


= A m *0 R A b _ , 


(A30) 


D (3) 


= 4 (0) *o^ , 


H (3) 


= Ai (0) *O fl ^ . 


(A31) 



Then the contributions of vertex loop diagrams (a,b,c) in Fig.l to the form factors and Hl,r are given as follows. 
• Diagram (a): 

F l ] = ; l 2 \-\ D(1) (- 1 + 4g 2 4 ) + m^orn B^C - m^m t C^ (C n - C 12 ) 

+m~-m t A^(C + C u - C 12 ) - m 2 t D^(C n - C 12 + C 21 - C 23 ) - m 2 w D^(C 22 - C 23 )} , (A32) 
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F R ] = -i I , \-\ E(1) (- 1 + 4C 24) + m^rn-- Co - m^m t F^ (C n - C 

AVtb I67H [2 x 3 x« *j 

+m.-m t H^(C + C n - C 12 ) - m 2 t E^(C n - C 12 + C 21 - C 23 ) - m 2 w E^(C 22 - C 23 )} , (A33) 
H ( L a) = ^[j^ {2m i oF^(C 11 - C 12 ) + 2m.-H^C 12 + 2m t E^(C n - C 12 + C 21 - C 23 )} , (A34) 
H R ] = 1^2 { 2m x° G(1) (Cii - C 12 ) + 2m--A«C 12 + 2m t D (1 \C 11 - C 12 + C 21 - C 23 )} , (A35) 

where the functions Co and C nm are 3-point Feynman integrals defined in Rcf. [34] with functional dependence 
Co , C nm (-p t , p w , m t ~ , m^o , m -- ) . 

• Diagram (b): 



W - ^! 1_/_R(2) 



L 2V tb 16 



(6) _ r_ L_ 

2F t6 16tt 2 



^{-5 (2) C 24 } , (A36) 
{-C (2) C 24 } , (A37) 



H ( ? = {m^o£»( 2 )(Co + C n ) + m t C^(C 12 - C n - C 21 + C 23 )} , (A38) 

^« ) = 2V tb lfar» { TO xS^ (2) (Co + Cn) + m t B^(C 12 - C n - C 21 + C 23 )} , (A39) 

with the Feynman integrals Co , C nm (—pt , j m x" > m t Q > m 6 (3 ) • 

• Diagram (c) : 

-^l^i ano " ^-^fl can ^ e obtained from the corresponding F^ R and if^jj with the replacements: 

4(1) ,4(3) , ^ B (3) 

C W _> C (3) , £>« -> r/ 3 > , - £( 3 ) , ff« - . (A40) 

For the self-energy loop diagrams in (d,e) of Fig. 1, we only present the corresponding forms of E^- (from chargino 
loops) and of E^o (from neutralino loops). The renormalization constants SZ^ b can be obtained in a straitforward 
manner by using Eqs.(6, 7). The results are given by 

• Diagram (d) : 

E**-(P) = ^^~ 2 {{\A\? $P L + \Atf p-PuWp,™--^) 

-m~- (A b *A\P L + A*?A b _P R )B (p, m-- , m Bfl )} , (A41) 

S£o(p) = {d4 0) | 2 ^ + |4 0) | 2 ^)Si(p,m^o,m t -J 

-n M (A!, 0) *4 0) P i + 4 0) *^ 0) Pfl)So(p,m i o ) m t -j} . (A42) 



Diagram (e) : 



-m--(A t *AVPL + ^A t _P fl ) J Bo(p,m i -,m f j} , (A43) 
- mi o(A m *Af ] P L + 4 (0) *4 6 _ (0) P fl ) J B (p,m s o,m 5(3 )} . (A44) 
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Comparing our results with those in [13], we found that our self-energy loop results essentially agree with those of 
[13], but our vertex loop results disagree with those of [13]. The analytic results of [13] seem to have some explicit 
errors, e.g., it seems to be impossible to cancel the UV divergence in Fl, terms like m t Co should not appear in Hl : r 
and G 2 and G3 should be interchanged. 



APPENDIX B: EXPLANATION OF THRESHOLD BEHAVIORS 



The two-point loop functions B' Q 1 (p,mi,m2) 



enter our results via SZj". In order to explain the threshold 



behavior in Figs. (2, 4, 6), we study the behavior of B' 1 (mo,mi,m2) near the threshold point m 
take B' Q as an example, the behavior of B[ is similar. 

For fixed mo and m 2 , B' (m , mi, m 2 ) as a function of mi can be expressed as 



mi + m 2 . We 



° Jo ™%x 2 - i m l 



x — x 



1 r 1 x—. 

m% JO (x — a) 



+ m\ — m 2 )x + ml 



dx 



dx 



where 



a = 



= < 



tUq + ml — m 2 
2mg 



(x-ayi-b^ 

tyIq JO {x — a) z 
1 rl x a; 2 7 

^ Jo { x - a y2 +V 2 ax 



for 
for 
for 



mi < mo - m 2 
mi = m — m 2 
mi > mo - m 2 



(Bl) 



2ml 



\Z\(m + mi + m 2 )(m + mi - m 2 )(m - mi + m 2 )(m - mi - m 2 )| . 



(B2) 
(B3) 



From Eq.(Bl), one may infer that, in regions below threshold (mi < mo — m 2 ) or above threshold (mi > mo — m 2 ), 
B' is an continuous function of mi. At the threshold point (mi = m — m 2 ), £> is divergent due to the fact < a < 1. 
(In some calculation programs, such as LoopTools [35] used in our calculations, an imaginary part is added to the 
propagators in the loops to avoid the divergence at the threshold point.) 

In the limit mi — > m — m 2 from the lower (i.e. mi < m — m 2 ), one has B' — * — 2 + (1 — 2a) In [(1 — a)/a] which 
is a finite value and corresponds to the top points in Figs. (2, 4, 6). But in the limit mi — > m — m 2 from the upper 
(i.e. mi > mo - m 2 ), one has £? — > 00 which corresponds to the bottom point in Figs. (2, 4, 6). Since these bottom 
points tend to take very large negative values, the perturbative calculation is not reliable at these points and a special 
treatment is needed [36]. 
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